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ABSTRACT 

Exoplanets in extremely close-in orbits are immersed in a local interplanetary medium 
(i.e., the stellar wind) much denser than the local conditions encountered around the 
solar system planets. The environment surrounding these exoplanets also differs in 
terms of dynamics (slower stellar winds, but higher Keplerian velocities) and am¬ 
bient magnetic fields (likely higher for host stars more active than the Sun). Here, 
we quantitatively investigate the nature of the interplanetary media surrounding the 
hot Jupiters HD46375b, HD73256b, HD102195b, HD130322b, HD179949b. We simu¬ 
late the three-dimensional winds of their host stars, in which we directly incorporate 
their observed surface magnetic fields. With that, we derive mass-loss rates (1.9 to 
8.0 x 10“ 13 M 0 yr -1 ) and the wind properties at the position of the hot-Jupiters’ or¬ 
bits (temperature, velocity, magnetic field intensity and pressure). We show that these 
exoplanets’ orbits are super-magnetosonic, indicating that bow shocks are formed sur¬ 
rounding these planets. Assuming planetary magnetic fields similar to Jupiter’s, we 
estimate planetary magnetospheric sizes of 4.1 to 5.6 planetary radii. We also derive 
the exoplanetary radio emission released in the dissipation of the stellar wind en¬ 
ergy. We find radio fluxes ranging from 0.02 to 0.13 mJy, which are challenging to 
be observed with present-day technology, but could be detectable with future higher 
sensitivity arrays (e.g., SKA). Radio emission from systems having closer hot-Jupiters, 
such as from r Boo b or HD 189733b, or from nearby planetary systems orbiting young 
stars, are likely to have higher radio fluxes, presenting better prospects for detecting 
exoplanetary radio emission. 

Key words: MHD - methods: numerical - stars: magnetic fields - stars: winds, 
outflows - planetary systems 


1 INTRODUCTION 

The interplanetary medium that surrounds exoplanets is 
filled by stellar wind particles and the embedded stellar 
magnetic field. Due to detection biases, the large major¬ 
ity of exoplanets found so far are orbiting cool stars at the 
main-sequence phase. Although the winds of these stars have 


proven quite challenging to observe (Mullan et al. 


1992 


Wargelin V Drake [2001 Wood et al. 2005), the interac¬ 


tion between exoplanets and their surrounding medium p.e., 
the host star’s wind) may give rise to observable signatures, 
such as planetary radio emission (Za rka|2007 ), enhancement 
of stellar activity (|Cuntz, Saar k, Mu sielak 2000; Shkolnik, 


Walker &; Bohlender|2003f Shkolnik et al.|2005 ), bow-shock 

formation (Vidotto, Jardine & Helling 2010 Llama et al. 


E-mail: Aline.Vidotto@unige.ch 


2013 Bisikalo et al.||2013 ), charge-exchange between stellar 
wind protons and planetary neutral hydrogen (Holmstrom 


et al.|2008| |Ekenback et al. 2010; Bo urrier k, Lecavelier des 


Etangs 2013; Kislyakova et al.|2014 ) and formation of comet¬ 


like tail structures ( Mura et al.|2011 Rappap ort et al.|2 012 

Budaj 2013| Villarreal D’Angelo et al. 2014), all of which 
can provide invaluable insights into the system, such as the 
intensity of the planetary magnetic field, velocity and tem¬ 
perature of the local stellar wind, etc. 

By studying stellar winds, we are able to make quantita¬ 
tive predictions about the interplanetary medium. A signif¬ 
icant improvement on our understanding of the interaction 
between a planet and the wind of its host star has been 
achieved in the past decade. Traditionally, these works have 
been based on simplified treatments of the winds (e.g. ip, I 
Kopp V Hu 2004| [Preusse et al.||2005| |GrieBmeier et al. 


2005 ; |Stevens |2005| Lovelace, Romanova &; Barnard||2008| 
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2 A. A. Vidotto et al. 


Vidotto, Jardine fe Helling||2011a| |See et aL]|2014| |Stru-| 
garek et al. |2014) . For example, simplified wind approaches 
might assume an isothermal wind structure, or that stars are 
non-rotating and/or non-magnetised bodies, among others. 
However, stellar winds are three-dimensional (3D) in na¬ 
ture, where complex interactions of a rotating, magnetised 
plasma take place. In view of that, more recently, new gener¬ 
ations of 3D, magnetohydrodynamical (MHD) models have 
started to be employed in the studies of interactions between 
stars/winds and their planets (e.g., | Vidotto et al.||2009a|b[ 


2010 [MTT1[MT2l[2014c||Cohen et al.|2009[|2014||Llama et al.| 


2013). 


The advantage of using simplified treatments of the 
wind is that these works rely on analytical and low¬ 
dimensional (ID, 2D) numerical studies, which are signif¬ 
icantly faster and do not demand extensive computational 
resources as 3D models do. Because of that, ID works can 
investigate a much wider range of stellar wind parameters 
(e.g. Cranmer & Saar|20TT) |See et al.|2014| ) than more com¬ 
plex, computationally-expensive 3D models can. The disad¬ 
vantage, on the other hand, is that the simplified models 
can not capture the 3D structure of stellar winds. Com¬ 
bined with modern techniques to reconstruct stellar mag¬ 
netic fields, some 3D models are able to provide a more 
realistic account of the stellar magnetic held topology em¬ 
bedded in the wind, recognised to be essential to interpret 


and predict signatures of star-planet interactions (|McIvor, 

Jardine & Holzwarth 

2006; [Fares et al. 2010; Vidotto, Jar- 

dine V Helling|]2011b 

Lanza 2012| Llama et al.|2013). 


1.1 Interactions with magnetised planets 

As the wind outflows from the star, it interacts with any 
planet encountered on its way. If these planets are magne¬ 
tised, their magnetic fields can act as shields, which prevent 
stellar wind particles from reaching all the way down to the 
surface or atmosphere of these objects (e.g. |Khodachenko 


et al.||20Q7| |Lammer et al.||2007| |Driscoll &; Bercovici||2013| 

Kisl yakova et al.|2014 ). This is in particular the case of the 

Earth and, more generally, of planets with dipolar held con¬ 
figurations. For these objects, the solar and stellar winds 
are dehected around the magnetospheric cavity, potentially 
helping the planet to retain its atmosphere. However, at¬ 
mospheric escape can still occur at high magnetic latitudes 


through polar hows, as is the case of the Earth (e.g., Seki 
et al.|[200T Moore & Horwitz 2007) and predicted for exo¬ 


planets (Owen &; Adams|2014 


see also Section 5.2). Part of 


this planetary outflow can return from the magnetosphere 
back into atmospheric regions of low-magnetic latitudes, re¬ 
ducing the total net loss rate of atmospheric escape, as sug¬ 


gested for the Earth scenario (Seki et al. 2001). The de¬ 


tailed process of atmospheric dynamics and escape is cer¬ 
tainly complex and not examined here. 

In the present work, only magnetised exoplanets are 
considered. This means that the cross-section of the ‘obsta¬ 
cle’ is not that of the planet itself, but rather takes into 
account the magnetospheric size of the planet. The magne¬ 
tospheric size of the planet depends both on the the charac¬ 
teristics of the local environment surrounding the planet (in¬ 
terplanetary density, velocity, magnetic field, temperature) 
and on its own magnetic held. On the theoretical side, some 
models suggest that the strength of the planetary magnetic 


held is dependent on the rotation rate of the planet (Farrell, 
Desch V Zarka 1999). In this situation, close-in planets that 


are tidally locked could have a reduced magnetic moment 


(Griehmeier et al. 2004). Other models advocate that the 


planetary magnetic held is related to the energy hux coming 
from the planetary core and does not depend on the rotation 
rate of the planet (Christensen, Holzwarth &; Reiners|2 009). 
Recent studies indicate that the planetary held strength is 
independent of rotation rate, which instead plays a role in 


the geometry of the generated magnetic held (Zuluaga V 
Cuartas||2012 ). 


Although planetary magnetism has been observed in 
several solar system planets, such as in Earth and the gi¬ 
ant planets, the presence of exoplanetary magnetic helds 


are much more elusive. Vidotto, Jardine V Helling (2010) 


suggested that the close-in giant planet WASP-12b hosts a 
bow-shock that surrounds its magnetosphere at a distance 
of about 4-5 planetary radii. Their suggestion was mo¬ 
tivated by transit observations of the close-in giant planet 


WASP-12b by Fossati et al. ([2010), who, based on space- 


borne spectroscopic observations in the near-UV, showed 
that the transit lightcurve of WASP-12b presents both an 
early ingress when compared to its optical transit, as well 


as excess absorption during the transit (see also Haswell 


et al. 2012). Vidotto, Jardine V Helling (2010) attributed 


this signature to an absorption of the material in the bow 
shock (see also Llama et al.||20lT ). If confirmed, this tech¬ 
nique should provide a useful tool for determining plane¬ 
tary magnetic field intensities for hot-Jupiter transiting sys¬ 
tems. In the case of WASP-12b, | Vidotto, Jardine V Helling| 
(2010) derived an upper limit of 24 G for the planetary field. 
Vidotto, Jardine V Helling (2011a) later proposed other 


targets with good prospects to hosting observable early- 
ingresses. Unfortunately, the near-UV (254 — 258nm) early- 
ingress signature of WASP-12b observed with (expensive) 
space-based spectroscopic observations ( |Fossati et al.||2010| 


Haswell et al. 2012) does not seem to be observable with 


ground-based, broad-band photometry in the wavelength 
range ~ 340 — 540nm (Ha swell et al.| |2012), and neither in 
the range of 303 — 417nm (Turner, private comm.; for other 
transiting exoplanets see | Turner et al.|2~0T3 ; Pearso n, Turner| 
V Sagan]|2014 ). Observations from Fossati et al.| ( |2010| ) in- 


dicate that the material surrounding WASP-12b absorbs at 
certain resonance lines in the near-UV (in particular in Mgll 
lines). The lack of absorption from broad-band photometric 
observations of WASP-12b possibly indicates that either the 
material is not absorbing at the observed photometric wave¬ 
lengths (~ 303 — 540nm), or that the absorption occurs only 
at some specific spectral lines, but gets diluted over the much 
wider spectral region. 

Another hint that close-in planets may also harbour in¬ 
trinsic magnetic fields, similar to the Earth and the giant 
planets of the Solar System, was found by Shkolnik, Walker 
V Bohlender (2003); Shkolnik et al. ( 2005| ~ 2008), who ob¬ 


served modulations of chromospheric spectral lines in phase 
with orbital periods on a few systems. Such modulations 
were interpreted as induced activity on the stellar surface 
due to magnetic interactions between star and planet. [Shkol-| 


nik et al. (2008) showed that there exists a correlation be¬ 


tween the night-to-night stellar activity variation with the 
ratio between the planetary mass to orbital period, used as 
a proxy for the magnetic moment of a tidally-locked planet. 
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On the environment surrounding close-in exoplanets 3 


Although unfortunately this correlation does not provide the 
intensity of the planetary magnetic field, it offers a way to 
measure the relative field strength among the different ex¬ 
oplanets in their sample. Therefore, once magnetism is as¬ 
sessed for one of their targets (by a different method), the 
magnetic field strength of their remaining targets could be 
derived. 

These two suggestions (early ingress and activity en¬ 
hancement), however, cannot be used as conclusive evidence 
of the presence of planetary magnetic fields, as alterna¬ 
tive, non-magnetic explanations for the observations exist 


(Preusse et al. 2006 

> Lai, Helling V van den Heuvel 2010 

Bisikalo et al. 2013| 

Vidotto et al. 2014a). A conclusive way 


to probe the presence of exoplanetary magnetic fields could 
be achieved by the detection of radio emission from the 
planet. The stellar wind that impacts on the planet pro¬ 
duces energetic particles that are captured by the planet’s 
magnetic field, emitting cyclotron radiation at radio wave¬ 
lengths. This emission depends on the planet’s magnetic 
field intensity and on the stellar wind power: it implies that 
the stronger is the stellar wind, the more luminous is the 
planet. As such radio emission is observed in the Solar Sys¬ 
tem (Zarka 2007), there are expectations that close-in exo¬ 


planets will exhibit comparable radiation (see |Nichols||2012 


for the case of planets that are not necessarily close-in). In 
particular, hot-Jupiters are expected to be much more lu¬ 
minous than the most luminous planet in our solar System, 
Desch &; Zarka|1999 ; GrieBmeier et al. 


Jupiter (e.g.. 


Farrell, 


2005 Zarka 2007; Jardine &; Cameron 2008; Vidotto et al. 


2010). This is because hot-Jupiters are located much closer 
to their stars, interacting with portions of the host-star’s 
wind that has larger kinetic and magnetic energies available 
to power planetary radio emission. So far, radio signatures of 
close-in exoplanets have not yet been det ected (e.g.|Bastian, 


Dulk &; Leblanc||2000| |Lazio et al.||2004| |Smith et al.| |2009 [ 

Halli nan et al.||2013| and one possible reason for that may 

be due to the lack of instrumental sensitivity in the appro¬ 


priate frequency range of the observations (Bastian, Dulk & 
Leblanc[?0 00 ). This picture, however, might be changing, as 
possible hints of exoplanetary radio emission have recently 


been reported (Lecavelier des Etangs et al. 2013 
et al.|[2014 ). 


Sirothia 


The theoretical estimates of the radio flux emitted by 
extrasolar planets carry a large uncertainty due to the fact 
that the stellar wind properties are poorly constrained. In 
this work, we model the 3-D structure of the stellar wind 
of a sample of 5 planet-hosting stars, whose 3-D winds have 
not yet been studied to date. We investigate the nature of 
the interplanetary media of these exoplanetary systems and 
how different they are from the environment surrounding 
our own Solar System planets. The stars used in this study, 
described in Section [ 2 ] have had their surface magnetic field 
recently reconstructed by means of tomographic techniques 
(Fares et al. 2012[ |2013| . These surface fields are used as 
boundary conditions for our data-driven simulations of stel¬ 
lar winds. Our model is described in Section |3] The derived 
global characteristics of the stellar winds are presented in 
Section [4] and the properties of the local environment sur¬ 
rounding the exoplanets in our sample are described in Sec¬ 
tion [5] We then use these computed quantities to calculate 
the strengths of the interactions between the stellar wind 
and the planetary system, making it possible to quantita¬ 


tively predict planetary radio emission and bow shock for¬ 
mation. Our discussion is shown in Section [ 5 ] and summary 
and conclusions are presented in Section [ 7 ] 


2 THE SAMPLE OF STARS 

The stars considered in this study consist of five solar- 
type stars of spectral types F 8 to Kl, namely: HD 46375, 
HD 73256, HD 102195, HD 130322 and HD 179949. All these 
stars host a gaseous planet at very close orbit (i.e., a hot- 
Jupiter). Table [l] presents a summary of the observationally- 
derived characteristics of the host stars and also of their hot- 
Jupiters (planet ‘b’). The large-scale surface magnetic field 
maps of the planet hosts have been reconstructed by | Fares] 
et al. ( 2012 , 20131 from a series of circular polarisation spec¬ 
tra (acquired at CFHT/ESPaDOnS and TBL/NARVAL) 
using the Zeeman-Doppler Imaging (ZDI) technique (e.g., 
Donati V Brown 199T| |Donati et al.] 2006). Figure [l] presents 
the radial component of the reconstructed surface field of 
these stars. Our targets present surface magnetic fields 
with a variety of topologies and intensities. For instance, 
HD 46375 presents a magnetic field that is mostly dipolar, 
whose axis is slightly tilted with respect to the rotation axis. 
HD 73256, on the other hand, has a magnetic field topology 
that is less axisymmetric. 


3 STELLAR WIND MODEL 

The stellar wind model we use here is identical to the 


one presented in Vidotto et al. (2014c). We use the 3D 


MHD numerical code BATS-R-US ( Powell et al.|1999 Toth 


et al.|2012 ) to simulate the stellar winds. BATS-R-US solves 
the set of ideal MHD equations for the mass density p, 
the plasma velocity u = {u r ,U0,U(p}, the magnetic field 
B = {B r , Be, Bcp}, and the gas pressure P : 


! + V-( P u) = 0, 


d(pu 

dt 


+ v. 


f T, 

V T BB 

r + ( p+ L 

) / _ ~ar_ 


de 

at 


ao 

—+ V-(uB-Bu) = 0, 


„ B 2 \ (u • B) B 


where 


e= P — + 

2 7 -I 


P B 
+ 8 ^' 


(1) 


= PS, ( 2 ) 


(3) 


= pg • u, (4) 


(5) 


We assume the wind is poly tropic, in which P oc p 1 and 7 is 
the polytropic index. To derive the temperature, we consider 
an ideal gas, so P — nksT, where ks is the Boltzmann 
constant, T is the temperature, n = p/(pm p ) is the particle 
number density of the stellar wind, pm p is the mean mass of 
the particle. In this work, we adopt 7 = 1.1, similar to the 


effective adiabatic index measured in the solar wind (Van 


Doorsselaere et al. 2011), and p = 0.5, for a fully ionised 
hydrogen plasma. 

At the initial state of the simulations, we assume that 
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4 A. A. Vidotto et al. 



Figure 1 . Mollweide projection of the radial component of the stellar surface magnetic field (from [Fares et a l. 2 012| 2 013| . These 
observed magnetic maps are included as boundary conditions in our data-driven simulations. 


Table 1. Observationally derived characteristics of the exoplanets and planet-host stars of our sample. The columns are: the host-star 
name, spectral type, mass (M*), radius (P*), effective temperature (T e ff), rotation period (P ro t), Rossby numbers (Ro), distance (d), 
inclination between the stellar rotation axis and the line-of-sight (z) estimated from ZDI, unsigned surface magnetic flux (4>o), date of 
the spectropolarimetric observations, projected mass of the planet c b’ (Mpsinz) and semi-major axis of the planetary orbit (a). For 
uncertainties in the quantities below, we refer the reader to the following literature. All the values listed below were compiled b y | Fares | 
|et al.| |2013| ), except for d , whose references are listed in the footnote of the table, and Ro, which was derived by |Vidotto et al.| ( |2014b| > 
using the models of |Landin, Mendes fe Vaz] ( |2010| >. 


Star 

ID 

Spectral 

type 

M* 

(Mq) 

R * 
(Re) 

T e ff 

(K) 

■Prot 

(d) 

Ro 

d 

(pc) 

z 

(deg) 

77 

o 

T-1 

Date 

M p sin z 
(M Jup ) 

a 

(«*) 

HD 46375 

K1IV 

0.97 

0.86 

5290 

42 

2.340 

33.4 1 

45 

0.85 

2008 Jan 

0.2272 

10.0 

HD 73256 

G8 

1.05 

0.89 

5636 

14 

0.962 

36.5 2 

75 

2.1 

2008 Jan 

1.869 

9.0 

HD 102195 

K0V 

0.87 

0.82 

5290 

12.3 

0.473 

29.0 3 

50 

2.1 

2008 Jan 

0.453 

12.6 

HD 130322 

K0V 

0.79 

0.83 

5330 

26.1 

0.782 

30.0 4 

80 

0.74 

2008 Jan 

1.043 

23.2 

HD 179949 

F8V 

1.21 

1.19 

6168 

7.6 

> 1.726 

27.0 5 

60 

1.3 

2007 Jun 

0.902 

7.9 


^Marcy, Butler &; Vogt| ( |2000| ; ^Udry et al.||2003| >; ^Ge et ah| ( |2006| >; ^Udry et al.||2000| ; ^Holmberg, Nordstrom &; Andersen| ( |2007| 


the wind is thermally driven (Parker 1958|. The stellar rota¬ 
tion period P ro t, M* and P* are given in Table[l] At the base 
of the corona (r m P*), we adopt a wind coronal tempera¬ 
ture To = 2 x 10 6 K and wind number density no = 10 9 cm -3 
(Section [6T] discusses the choices of no and To and how they 
affect our results). With this numerical setting, the initial 
solution for the density, pressure (or temperature) and wind 
velocity profiles are fully specified. The radial component 
of the magnetic field B r anchored at the base of the wind, 
is reconstructed from observations (Fig. [l]). The other two 
components of the surface field are assumed to be poten¬ 
tial (V x B = 0), as it has been shown that stellar winds 
are largely unaffected by the non-potential part of the ob¬ 
served surface field ( |Jardine et al.|2013 ). At the initial state, 
we assume that the field considered in the simulation box 
is potential up to a radial distance r = rss (known as the 
source surface) and, beyond that, the magnetic field lines 
are considered to be open and purely radial. As the simu¬ 
lation evolves in time, the wind particles interact with the 
magnetic field lines (and vice-versa), removing the field from 
its initial potential state. For all the cases studied here, we 
take rss = 4 P*, but we note that different values of rss pro¬ 
duce similar final steady-state solutions for the simulations 
( Vidotto et al.|20lT 2014c). 


Once set at the initial state of the simulation, the val¬ 


ues of the observed B r are held fixed at the base of the 
wind throughout the simulation run, as are the coronal base 
density and thermal pressure. A zero radial gradient is set 
to the remaining components of B and u = 0 in the frame 
corotating with the star. The outer boundaries at the edges 
of the grid have outflow conditions. The rotation axis of the 
star is aligned with the z-axis, and the star is assumed to 
rotate as a solid body. Our grid is Cartesian and the star is 
placed at the origin of the grid, which extends in x, y, and z 
from —20 to 20 P*, except for HD 102195, whose simulation 
box extends from —24 to 24 P*, as to extend out to the 
orbit of the planet. BATS-R-US uses block adaptive mesh 
refinement. The finest resolved cells are located close to the 
star (for r < 2 P*), where the linear size of the cubic cell 
is 0.0097 P* (or 0.012 P* for the simulation of HD 102195). 
The coarsest cell has a linear size of 0.31 P* (or 0.37 P* for 
HD 102195) and is located at the outer edges of the grid. 
The total number of cells in our simulations is around 40 
million. As the simulations evolve in time, both the wind 
and magnetic field lines are allowed to interact with each 
other. The resultant solution, obtained self-consistently, is 
found when the system reaches steady state in the reference 
frame corotating with the star. 
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4 DERIVED PROPERTIES OF THE STELLAR 
WINDS 

Table [2] presents the properties of the stellar winds obtained 
in our simulations. The unsigned observed surface magnetic 
flux is 

3>o = j \B r (R*,6,<p)\dS* (6) 

and the unsigned open magnetic flux is 

^open — f \B r (r,e,v)\dS sph . ( 7 ) 

'^'S'sph 

The surface flux (Table [lj is integrated over the surface of 
the star S'* and the open flux (Table [ 2 ]) over a spherical 
surface S sp h at a distance r from the star, where all the 
magnetic held lines are open. The mass-loss rate M of the 
stellar wind, which outflows along open magnetic held lines, 
can be calculated as the hux of mass integrated across S sp h 

M — (D pU r dS S phi ( 8 ) 


where M is a constant of the wind. Similarly, the angular 
momentum loss rates can be calculated as the angular mo¬ 
mentum hux across S sp h 


J = 


'uj -B (p 5 ^ 

47T 


+ mup>pu r 


d*Ssph 


( 9 ) 


( Mestel &; Selley|[l970| Mestel|p~999| Vidotto et al. 2014c), 


where vj — (x z +y z ) l/z is the cylindrical radius. In our simu¬ 
lations, we hnd that M ranges from ~ 2 to 8 x 10“ 13 M© yr -1 
and J between ~ 0.14 and 2.4 x 10 31 erg for the stars in 
our sample. The open hux ranges from 26% to 69% of the 
large-scale unsigned surface hux. These values are within 
the range (4.4 — 8.4) x 10 22 Mx. For the solar wind, Wang, 
Sheeley & Rouillard ( |2006 ) obtained magnetic held values 
at the orbit of the Earth in the range between 0.01 and 
0.05 mG, or in terms of open magnetic huxes, in the range 
(2.8— 14) x 10 22 Mx, depending on the phase of the solar ac¬ 
tivity cycle. Although the range of open huxes calculated for 
the simulations presented here fall within the values of the 
solar wind, we show in Section [5] that the values of the in¬ 
terplanetary magnetic held at the orbits of the hot-Jupiters 
are more than 100 times larger than the interplanetary mag¬ 
netic held at the Earth’s orbit (compare 0.01 to 0.05 mG to 
the values presented in Table j3b. 

The left panels in Figure [2] show the hnal conhguration 
of the magnetic held lines obtained through self-consistent 
interaction between magnetic and wind forces after the sim¬ 
ulations reached steady state. Although we assume the mag¬ 
netic held is current-free in the initial state of our sim¬ 
ulations, this conhguration is deformed when the interac¬ 
tion of the wind particles with the magnetic held lines (and 
vice-versa) takes place (currents are created in the system). 
The right panels of Figure [2] show the Alfven surface Sa 
of each simulation. This surface is dehned as the location 
where the wind velocity reaches the local Alfven velocity 
(va — R(47rp) _1//2 ). Inside Sa, where the magnetic forces 
dominate over the wind inertia, the stellar wind particles 
are forced to follow the magnetic held lines. Beyond Sa, 
the wind inertia dominates over the magnetic forces and, 
as a consequence, the magnetic held lines are dragged by 


the stellar wind. In models of stellar winds, the Alfven sur¬ 
face has an important property for the characterisation of 
angular momentum losses, as it dehnes the lever arm of the 


torque that the wind exerts on the star (e.g., Weber Sz Davis 
1967). Its location is also relevant in the studies of magnetic 


interactions with planets (e.g., Strugarek et al. 


hen et al. 2014). As shown in Vidotto et al. (2014c), the 


2014 


Co- 


Alfven surfaces of the objects investigated here have irreg¬ 
ular, asymmetric shapes as a consequence of the irregular 
distribution of the observed magnetic held. To illustrate the 
difference in sizes of these surfaces, we show in the right pan¬ 
els of Figure [5] the scales of the images plotted (red lines). 
We hnd that the average radius of the Alfven surfaces range 
between 2.8 R* (for HD 179949) to 6.4 R* (for HD 102195). 

In order to provide constraints for analytical methods of 
extrapolation of magnetic held lines, we also compute here 
the MHD equivalent of the source surface. In particular, the 
potential held source surface (PFSS) method has proven to 
be a fast and simple way to extrapolate surface magnetic 
helds into the stellar coronal region (|Jardine et al.|1999 Jar- 


dine, Collier Cameron &; Donati|2002 Vidotto et al. 2013). It 
is also used here as the initial conditions for our simulations. 
However, the PFSS method has an unconstrained parame¬ 
ter: the radius rss of the source surface, beyond which the 
magnetic held lines are assumed open and purely radial, as a 
way to mimic the effects of a stellar wind. Because of stellar 
rotation and magnetic held stresses, in the MHD solutions, 
the surface where all magnetic held lines are purely radial 
does not exist - even in the region of open held lines, there is 
always Be and, especially, B^ components that are non-null. 
Therefore, we dehne here an “effective radius of the source 
surface” r|f as the radius of the spherical surface where 
97 percent of the average magnetic held is contained in the 


radial component (i.e., (|R r |)/(|R|) = 0.97, based on Riley 
et al.| 2006). For some of the stars in our sample (HD 73256 
and HD 179949), the ratio (\B r \)/(\B\) does not reach the 
97-percent level and in such cases, we take f to be the po¬ 
sition where (|R r |)/(|R|) is maximum. Table [ 2 ] shows that 
r|f is in the range between 1.9 R* and 5.6 R*, indicating a 
compact region of closed held lines. We note that this size 
is similar to the usual adopted size of 2.5 R© from PFSS 
methods of the solar coronal magnetic held and also similar 
to the values obtained in other MHD simulations of winds 
( Riley et al.||2006 Vidotto et al.|20lT 2014c). 


5 CHARACTERISING THE LOCAL 
ENVIRONMENT SURROUNDING 
HOT-JUPITERS AND RESULTANT 
INTERACTIONS 

All the stars in our sample host giant planets orbiting at 
close distances. Mercury, the closest planet to our Sun, has 
a semimajor orbital axis of about 0.39 au, or equivalently, of 
about 83 R©. The hot-Jupiters in our sample have consider¬ 
ably closer orbits, with semimajor axes of about 9 to 23 R* 
(i.e., about 9 to 4 times closer than Mercury). As a conse¬ 
quence, the hot-Jupiters in our sample interact with much 
denser winds that have larger ram pressures than those typ¬ 
ically found around the planets in the solar System. In addi¬ 
tion, because the hot-Jupiters are located much closer to the 
star, the large-scale magnetic field at the orbit of these plan- 


© 2014 RAS, MNRAS 000,|THT5l 

















































6 A. A. Vidotto et al. 


Table 2. Characteristics of the stellar winds. The columns are: the star name, the stellar wind mass-loss rate (M), angular momentum- 
loss rate (J), unsigned open magnetic flux ($ 0 pen), average radii of the Alfven surfaces ((r^)), its minimum and maximum value (r™ m 
and r™ ax ) and the effective radius of the source surface derived from the MHD models (rjf). In our simulations, M, J and <I>open are 
conserved within 0.03%, 3% and 4%, respectively. 


Star 

ID 

M 

(10- 13 Mq yr- 1 ) 

j 

(10 31 erg) 

^open 

(*o) 

{ta) 

(ft*) 

r„min max] 

l r A » r A J 

(ft*) 

r eff 

r SS 

(ft*) 

HD 46375 

1.9 

0.14 

0.52 

5.1 

[ 3 . 0 , 6 . 1 ] 

2.7 

HD 73256 

2.1 

2.3 

0.26 

6.2 

[1.8,8.1] 

5.6 

HD 102195 

3.2 

2.0 

0.41 

6.4 

[2.3, 7.5] 

5.6 

HD 130322 

5.8 

0.36 

0.69 

3.5 

[1.6,4.2] 

1.9 

HD 179949 

8.0 

2.4 

0.34 

2.8 

[1.0,3.7] 

3.0 


ets has also a larger strength compared to the interplanetary 
magnetic field strength of our solar System planets. 


The orbital planes of the planets considered in this work 
are not known. Here, we assume their orbits lie in the equa¬ 
torial plane of the star. This seems to be a reasonable hy¬ 
pothesis for our targets (cf. Table [l|, as planets orbiting 
stars cooler than 6200 K have been observed to have small 
(projected) obliquities ( |Winn et al.||201 0). Figure [3] shows 
the total pressure ptot (be., the sum of thermal, magnetic 
and ram pressures) experienced by a planet as it orbits at 
the equatorial plane of the stars. Note that the ram pres¬ 
sure term must take into account the relative motion of the 
planet through the interplanetary medium. Here, we assume 
prograde motion of the planetary orbit relative to the stel¬ 
lar rotation. The white circles indicate the orbital radii of 
each hot-Jupiter, taken here to be circular (note that for 
the systems investigated here the eccentricities are rather 
small, < 0.06). The colour-bar is the same for the five im¬ 
ages, illustrating that the total pressure varies from planet 
to planet. The last panel in Figure [3] shows the total local 
pressure at the planetary orbits as a function of subplane¬ 
tary longitude (see also Table [3]). For the cases studied here, 
at these orbital distances, the dominant term in the total 
pressure is the ram pressure of the relative motion of the 
planet through the wind. The values of the local total pres¬ 
sure are within (0.58 — 4.1) x 10 -4 dyn cm -2 , which are 
about 4 orders of magnitude larger than the ram pressure 
of the solar wind at the Earth’s orbit (1.8 x 10 -8 dyn cm -2 , 
See et ah] 2014). We also note that there is some variability 
in the local total pressure, showing that the planets interact 
with the varying environment of the star along their orbits. 
In the case of HD 73256, the amplitude of this variability 
is the highest among the cases studied here and is due to 
the peak (which is a factor of 1.9 above the average value 
of ptot of HD 73256) at ~ 220 deg. This peak is caused 
by a fast wind stream, associated to the magnetic feature 
seen in the surface magnetograms at longitude ~ 225 deg 
(Fig. [l}. A similar feature appears in Figures [4] and [5] that 
we present later. Variability on larger timescales due to in¬ 
trinsic variations of the stellar magnetic field can also alter 


the environment surrounding planets (Vidotto, Jardine V 


Helling 2011b Vidotto et al.||20l2 Llama et al.||2013 ), but 

it is not considered in the present work. 


5.1 Exoplanetary bow shocks: sizes and 
orientations 

If a planet is magnetised, its magnetic field can act as shield 
for the stellar wind, deflecting the wind particles and po¬ 
tentially preventing the wind from reaching down to the 
planetary atmosphere. A way to estimate the size of this 
stand-off distance is by pressure balance between the local 
total pressure of the interplanetary medium (i.e., the stellar 
wind) and the planet total pressure. Thus, at the interaction 
zone, we have 


Bi 


Ptot - 


87r 


( 10 ) 


where B P ^ M is the planetary magnetic field intensity at a 
distance vm from the planet centre. Eq. (10) neglects the 


planetary thermal pressure component on the right side. Be¬ 
cause of the exponential decay of planetary densities, at the 
height of a few planetary radii, the thermal pressure is usu¬ 
ally negligible compared to the planetary magnetic pressure. 
If we assume the planetary magnetic field is dipolar, we have 
that B p , rM = B p , ec[ (R p /rM ) 3 , where R p is the planetary ra¬ 
dius and B Pj e q its surface magnetic field at the equator (half 
the value of the intensity at the magnetic pole). For a plan¬ 
etary dipolar axis aligned with the rotation axis of the star, 
the magnetospheric size of the planet is given by 


tm 

Rv 


Bi 


87rp t ot 


1/6 


(ii) 


In the absence of observational constraints, we assume the 
hot-Jupiters studied here to host magnetic fields similar to 
Jupiter’s. Figure [4^t shows the magnetospheric sizes of these 
hot-Jupiters assuming q = 7 G (i.e., half of the maxi¬ 
mum observed field of Jupiter of ~ 14 G, Sm ith et al.|1 975; 
Bagenal 1992 ). The average estimated magnetospheric sizes 
range from about (t*m) = 4.2 R p for HD 179949b to ( t*m ) = 
5.6 R p for HD 130322b (see Table[3]). Variations in vm along 
the planetary orbit are roughly ~ 10%. This variation oc¬ 
curs because, as the planet goes along its orbit and probes 
regions with different ptot, its magnetospheric size reacts ac¬ 
cordingly, becoming smaller when the external ptot is larger 
and vice-versa. 

Over-plotted to Figure [3] are the contours at the equa¬ 
torial plane of the Alfven surface (red lines) and the mag- 
netosonic surface (black lines) of the stellar wind. In all the 
cases studied here, the planets orbit at regions of super fast 
magnetosonic velocities of the stellar wind. One exception 
is the case of HD 73256b, in which a small part of its or- 
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Figure 2. Left: The final configuration of the magnetic field lines after the wind solution has relaxed in the grid. Over-plotted at the 
surface of the star is the observationally reconstructed stellar magnetic field ( |Fares et a l. 2012] [2013]), used as boundary condition for 
the radial magnetic field. Right: The Alfven surfaces are shown in grey. Note their irregular, asymmetric shapes due to the irregular 
distribution of the observed field. The equatorial ( xy ) planes of the star, assumed to contain the orbits of the planet, are also shown, as 
are the intersections between the xy plane and the Alfven surface (thin black contour) and the orbital radius of the planet (thick blue 
contour). 


bit (white circle) lies within the fast magnetosonic surface 
of the wind. This does not necessarily mean that at these 
orbital positions a bow shock will not be formed surround¬ 
ing HD 73256b’s magnetosphere. Rather, it is the relative 


velocity of the planet orbiting through the stellar wind 

Au = u — uk<P, (12) 

where uk is the (purely azimuthal) Keplerian velocity of the 
planet, that should be compared to the fast magnetosonic 
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HD130322 


20 Rstar 


20 Rstar 


HD179949 


Figure 2 - continued 


Table 3. Derived characteristics of the hot-Jupiters and of their local environments. The columns are, respectively, the planet name, the 
averages of the local velocity of the wind in the reference frame of the planet, wind density, magnetic field strength, wind temperature, total 
pressure, planetary magnetospheric radius, shock angle, auroral oval opening angle and fractional area of the polar cap. These quantities 
were averaged over the subplanetary longitude. Values in brackets represent the minimum and maximum values of the averaged quantity. 


Planet 

ID 

(Au) 

(km s -1 ) 

<n) 

(10 5 cm -3 ) 

<|B|> 

(mG) 

<T) 

(10 6 K) 

(ptot ) 

/IQ- 4 dyn ) 
v cm z J 

( r l Vf) 

(R P ) 

($shock) 

(deg) 

<«o) 

(deg) 

(^auroral) 

(^planet) 

(0radio) 

(mjy) 

HD46375b 

234 

1.8 

8.8 

0.87 

l.i 

5.1 

52 

26.2 

0.10 

0.037 


[228, 242] 

[1.7, 2.0] 

[0.55,11] 

[0.86,0.91] 

[1.0,1.2] 

[5.0, 5.2] 

[50, 53] 

[26.0, 26.6] 

[0.10,0.11] 

[0.036, 0.043] 

HD73256b 

263 

2.0 

17 

1.1 

1.6 

4.8 

57 

27.1 

0.11 

0.045 


[217, 345] 

[1.6, 2.6] 

[2.6, 26] 

[0.91,1.6] 

[1.0, 2.7] 

[4.4, 5.2] 

[44,67] 

[26.1,28.5] 

[0.10,0.12] 

[0.027, 0.081] 

HD102195b 

288 

1.5 

14 

0.96 

1.3 

5.0 

65 

26.6 

0.11 

0.067 


[240, 338] 

[1.1, 2.0] 

[3.6,18] 

[0.87,1.2] 

[1.1, 1.6] 

[4.8, 5.1] 

[61,69] 

[26.2,27.1] 

[0.10,0.11] 

[0.054, 0.086] 

HD130322b 

322 

0.6 

2.3 

0.78 

0.62 

5.6 

74 

25.0 

0.09 

0.055 


[316, 334] 

[0.6, 0.7] 

[0.36,2.9] 

[0.77,0.79] 

[0.58, 0.69] 

[5.5, 5.7] 

[74,75] 

[24.8, 25.2] 

[0.09,0.10] 

[0.053, 0.061] 

HD179949b 

243 

5.9 

9.6 

0.97 

3.8 

4.2 

53 

29.3 

0.13 

0.112 


[225, 257] 

[5.5, 6.2] 

[1.4,15] 

[0.96,0.99] 

[3.1,4.1] 

[4.1,4.3] 

[51, 55] 

[28.9, 29.6] 

[0.12,0.13] 

[0.092,0.127] 


velocity of the local plasma Vf — (c 2 + v^) 1 ^ 2 , where c s 
is the local sound speed. Figure [4 ]d shows the fast magne- 
tosonic Mach number (A u/vf) calculated at the orbital radii 
of the hot-Jupiters, where we see that the relative planetary 
velocity is always super-fast magnetosonic (i.e., A u/vf > 1 ), 
indicating that the magnetosphere of these planets are sur¬ 
rounded by bow shocks. 


It has been proposed that these bow shocks might ab¬ 
sorb at specific wavelengths, generating asymmetric transit 
lightcurves (Vidotto, Jardine <V Helling||2010 ). This is par¬ 
ticularly relevant for the case of hot-Jupiters, in which the 
orientation of the bow shock is shifted towards the direction 
of planetary motion (as opposed to the bow shocks surround¬ 
ing the solar system planets, which are largely formed facing 
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Figure 3. Distribution of the total pressure ptot experience by a planet as it orbits at the equatorial plane of each star in our simulations. 
The black dashed and solid red contours are cuts of the fast magnetosonic and Alfven surfaces of the stellar wind, respectively, at the 
equatorial plane. The white lines indicate the orbital radii, taken to be circular, of the hot-Jupiters. The last panel shows the total local 
pressure at these orbits as a function of subplanetary longitude. 


the Sun). These ‘sideways’ bow-shocks present the best con¬ 


ditions for detection during planetary transits (Llama et al. 


2011 


2013). Although the hot-Jupiters investigated here are 


not transiting and do not have constrained orbital inclina¬ 
tions, there has been cases in the literature of non-transiting 
(but grazing) exoplanets whose extended atmospheres might 
undergo partial transit (Ehrenreich et al. 2012). Likewise, 


combined with the orbital inclinations, it is possible that 
bow shocks of non-transiting planets might be visible if they 
graze the stellar disc. We, here, do not model the 3D extent 


of bow shocks, as done in Llama et al. (2011 2013), but we 


can calculate the angle between the shock normal and the 
tangent of a circular orbit 


^shock — arctan 


\Uk — u< 


(13) 


(Vidotto, Jardine V Helling 2010). Along its orbital path, 


the planet probes regions of the wind with different veloc¬ 
ities, which implies that the orientation of the bow shock 
that forms surrounding planetary magnetospheres changes 
along the planetary orbit. This can be seen in Figure H 
and Table [3] where we present 0 s hock as a function of the 
subplanetary longitude. We present in Table [ 3 ] the average 
shock angle (0 s hock) of the bow shock of each of these hot- 
Jupiters, where we note that they range from about 52° to 
74°. 


5.2 Exoplanet ary auroral ovals: escape channels 
and radio emission 

As the planetary magnetosphere extent is reduced, the size 
of the ‘auroral oval’, which is the amount of planetary area 
with open magnetic field lines, increases. Along these open 
field lines, particles can be transported to/from the inter¬ 
planetary space, affecting, for instance, the amount of at¬ 


mospheric mass loss (Adams 2011). We estimate here the 


size of the auroral region of the planet as follows. Assum¬ 
ing the planet to have a dipolar magnetic field, aligned with 
the planetary orbital spin axis, the colatitude of the largest 
closed field line of the planet, which defines the boundary be¬ 
tween open- and closed-field line regions, can be estimated 


as ao = arcsin( J R p /rM) 1//2 (Siscoe V Chen 


1975 


Tarduno 


et al.|[201() ). This implies in a fractional area of the plane¬ 


tary surface that has open magnetic field lines 

-^polar cap \ 

— A -- = (1 — cos ao), 

-^■planet 


(14) 


( Vidotto et al.| 2013). Therefore, in addition to making tm 
smaller, a stronger external pressure of the stellar wind ex¬ 
poses a larger area of the polar cap of the planet. Table [ 3 ] 
shows the average, minimum and maximum angles of the 
auroral ovals (ao) and fraction of open area (H po i ar cap ) as 
calculated by Eq. ( fl4| . For the hot-Jupiters analysed here, 
(ao) ranges between 25° and 29°, and (A po iar cap) ranges be¬ 
tween 9% and 13%. For comparison, the size of the auroral 
oval in Saturn is ao — 10° - 20° ( | Clarke et al. 2005) and 
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longitude (°) 



longitude (°) 



longitude (°) 


Figure 4. As a magnetised planet orbits around its host star, 
it probes regions of the stellar wind with different properties. 
As a consequence, its magnetospheric size and shock orientation 
change. Upper panel: the magnetospheric stand-off distance for 
the hot-Jupiters studied here as a function of subplanetary longi¬ 
tude. Middle panel: the ratio between the relative velocity of the 
planet and the local fast magnetosonic velocity. Bottom panel: the 
angle formed between the shock normal and the tangent of a cir¬ 
cular orbit. The top and bottom figures assume the hot-Jupiters 
have a dipolar field of 7 G at their equator. 


at the Earth it is ao ~ 17° - 20° (Milan et al. 2009). Using 
Eq. (14), a rough estimate indicates that the open-field-line 
region covers ~ 1.5% - 6% of Saturn’s surface and ~ 4.5% - 
6% of Earth’s surface. This is a factor of ~ 2 smaller than the 
values we derive for the hot-Jupiters in our sample, but not 
as extreme as the cases of planets orbiting at the habitable 
zone of more active stars (Vidotto et al. 2013). 

Planetary radio emission takes place in a hollow cone of 
half-aperture angle given by the auroral oval co-latitude ao • 
It has been recognised that the radio emission of the Earth 
and the four giant planets of the solar system correlates to 
the local characteristics of the solar wind (e.g., Zarka [l998 ), 
an indication that radio emission is powered by the local 
solar wind. Analogously, it is expected that when exoplanets 
interact with the wind of their host stars, they would also 
be sources of radio emission. 

We use the results of our stellar wind simulations to 
calculate the kinetic power of the wind, at the orbital radii of 
the hot-Jupiters studied here. Our approach follows closely 


the one in Vidotto et al. (2012). The kinetic power Pk of the 


impacting wind on the planet is approximated as the ram 
pressure of the particles p(Au) 2 impacting on the planet, 
with effective cross-section 7rr|^, at a relative velocity Au 


Pk ~ p(Au) 3 7vrM. 


The radio flux can be written as 

Pradio 


0radio — 


VkPk 


d 2 uAf d?uAf 


(15) 


(16) 


where d is the distance to the system, uj — 2 x 2n(l — cos ao) 
is the solid angle of the hollow emission cone (defined by 
the auroral oval), and Af is the frequency of emission. In 
the last equality, we assumed a linear efficiency rjk in con¬ 
verting the power released from the dissipation of kinetic 
wind energy to radio emission (‘radiometric Bode’s law’). 
We adopt 77 k = 10 -5 , as derived from observations of the 
Solar System planets (Zarka [2007). Here, we assume that 
the emission bandwidth Af is approximately the cyclotron 
frequency ( Griefimeier et aL]|2007| : 


Af = e JP°Al — 2.8 


m e c 


Bp (c^o) 


1 G 


MHz, 


(17) 


where m e is the electron mass and c the speed of light. 
B p (a 0 ) is the planet’s magnetic field strength at colati¬ 
tude ao of the auroral ring. For a dipolar field, B p (a 0 ) = 
B p ; eq(1 + 3 cos a 0 ) 1/2 . 

To compute the radio flux (Eq.[l6]), we need to know the 
physical size of vm- This value, normalised to the planet’s 
radius, is given in Figure [4 ^l and we further assume plane¬ 
tary radii of 1.5i?j up for all the hot-Jupiters analysed in this 
work (note that they are non-transiting planets and there¬ 
fore do not have observationally-determined radii). Eq. (16) 
is the only place where the physical size of the exoplanet is 
required and different choices of R p influence the estimated 
radio flux as </> r adio oc r 2 M oc R 2 . 

Figure [5^i shows the radio flux computed using the re¬ 
sults of our wind simulations and Figure [5}) shows the cal¬ 
culated frequency of emission. We find that the predicted 
emission frequency occurs at ~ 36 MHz and the radio fluxes 
range between 0.02 and 0.13 mJy among all the cases stud¬ 
ied here (see also Table [3| . Values of radio fluxes such as 
these (including the peak values that occur at favourable 
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longitude (°) 


Figure 5. The predicted radio flux (Eq. |16| > computed using the 
results of our wind simulations (top) and associated frequency 
of emission (bottom) assuming the emission bandwidth is the 
cyclotron frequency (Eq. ED- These results assume a dipolar ex¬ 
oplanetary magnetic field, whose intensity is 7 G at the equator. 


phases) should be challenging to be observed with present- 
day technology, such as with LOFAR, whose sensitivity at 
20 to 40 MHz is > 30 to 3mJy, respectively, for a one-hour 
integration time ( |Grie fimeier, Zarka & Girard 2011). It is 
likely, however, that even these small radio fluxes will be 
detectable with future higher sensitivity arrays, such as the 
SKA-low array system. 

Among the systems studied here, HD 179949b has the 
highest estimated radio flux. This occurs for two reasons. 
First, this exoplanet has the closest orbital radius and, be¬ 
cause of that, pAu 3 is the largest among our sample; for the 
same reason, it also has the smallest vm (cf. Tables [l] and |3|. 
In spite of the smallest cross-section tv r 2 M , the large pAu 3 
term is more important in Eq. (15), which results in the 
largest stellar wind kinetic power impacting on the magne¬ 
tosphere of the exoplanets studied here. Second, the closest 
distance to the HD 179949 system also favours a larger radio 
flux (Eq. (16)). 

It is also worth comparing the emission calculated here 


and the values calculated for r Boo b and HD 18973310 Us¬ 
ing the same radio emission model presented here, |Vidotto| 


et al. (2012) estimated the radio flux of r Boo b at different 


epochs of the host star’s magnetic cycle. They found the ra¬ 
dio flux of r Boo b to be of the order of 0.5 - 0.9 mJy. We 
can also use the simulations presented in Llama et al. (2013) 
to compute the radio flux of HD 189733b. Assuming a plan¬ 
etary radius of R p = 1.15 Rj U p and a distance of 19.3 pc, 
we calculate the radio flux of HD 189733b to be on average 
0.47 mJy (peak at 0.98 mJy) for the case where the observed 
stellar magnetic map is derived from the 2007 June obser¬ 
vations and 0.23 mJy (peak at 0.51 mJy) for the 2008 July 
map (cf. |Fares et al. 2010J. 

The radio fluxes computed for r Boo b and HD 189733b 
are therefore considerably larger than the values computed 
for the exoplanets presented here, having better prospects 
for being detected. The reason why these two systems have 
higher radio fluxes is similar to the reasons discussed for 
the case of HD 179949b: a combination of closer orbital 
radii (6.8 R* and 8.6 R* for r Boo b and HD 189733b, 
respectively) and closer distances to the systems (15.6 and 
19.3 pc). It is also expected that exoplanets orbiting young 
stars (with denser stellar winds) are likely to produce higher 
radio fluxes ( Griefimeier et al.|2005| Vidotto et al.|2010 ), pre¬ 
senting also better prospects for detection of exoplanetary 
radio emission. 

Radio fluxes estimated for the 5 hot Jupiters stud¬ 
ied here have been estimated by other authors. For in¬ 
stance, Griefimeier, Zarka & Girard (2011) predicted ra¬ 


dio fluxes that are larger than the values predicted here 
by a factor of 500 - 2000 (compared to the case for their 
rotation-independent planetary magnetic field model). Al¬ 
though our radio emission model is similar to the one used 
Griefimeier, Zarka &; Girard (201l| (i.e., both our mod¬ 


els assume a ‘radiometric Bode’s law’, in which a fraction of 
the dissipation of the wind power is converted into planetary 
radio emission), we attribute the difference found between 
their work and the present one due to the different mod¬ 
els assumed for the stellar wind and stellar magnetic field, 
as well as for the assumed planetary magnetic field intensi¬ 
ties. For the stellar wind, [Griefimeier, Zarka fc Girard) s work 
assumes a spherically symmetric, isothermal wind model 
(Parker 1958). The velocity and density structures are scaled 
with respect to the age of the system, based on the age re¬ 
lations found by Newkirk (1980) and Wood et al. (2005). 


For the planetary magnetic field, |Griefimeier, Zarka fe Gn] 
|rard[ s work assumes either a case where the planetary dy¬ 
namo is independent (Reiners & Christensen 2010) or de¬ 
pendent (Griefime ier et al.||2004 ) on the planetary rotation. 


Griefimeier, Zarka & Girard (2011) showed that the inten¬ 


sity of the planetary magnetic field affects the frequency of 
the emission (as in our model) and that the radio flux has a 
strong dependence with the intensity of the planetary mag¬ 


netic field (contrary to our model). More recently, See et al. 


(2015) studied the variability of exoplanet ary radio emis¬ 
sion for a sample of planet host stars, which includes the 
objects studied in the present work. Similar to our model, 


1 Note that the stellar wind simulations presented in |Vidotto| 
|et al.| p012| ) and |Llama e t al. (2013) have the same assumptions 
as the ones shown in the present work. 
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their model incorporates the realistic large-scale geometry of 
the stellar magnetic field, but their radio emission model dif¬ 
fers from ours. Instead, their study was based on the model 


developed by Jardine k Cameron (2008), which computes 


the radio emission generated by energetic electrons released 
in the reconnection between stellar and exoplanetary mag¬ 
netic field lines, without assuming the a priori relation of 
the ‘radiometric Bode’s law’. Despite the differences in these 


models, the radio fluxes estimated by See et al. (2015) and 


by us are very similar, within a factor of 1 - 4, except for 
HD 130322, in which we estimate radio fluxes that are 100 
times larger than theirs. In addition to providing informa¬ 
tion on exoplanet ary’s magnetic field, detection of exoplan¬ 
etary radio emission would clearly provide invaluable con¬ 
straints to stellar wind models as well. 


6 DISCUSSION 

6.1 Limitations of the models 

The stellar wind models presented in this paper use as in¬ 
put the observationally reconstructed stellar magnetic field 
and are, therefore, more realistic (and provide an advance) 
compared to models that are non-magnetised or that as¬ 
sume simplified stellar magnetic field topologies. In spite 
of that, our wind models share the limitations of global, 
polytropic wind models. In particular, these types of models 
have three parameters that are poorly constrained by ob¬ 
servations, namely, the wind base density and temperature 
and the temperature profile (i.e., the profile of energy de¬ 
position through the parameter 7 ). In this work, we have 
chosen to set all these three parameters to be the same for 
all the stars in our sample. On the other hand, parameters 
such as the stellar mass, radius, rotation period and mag¬ 
netic field differ for each object and are constrained to values 
observationally-derived for each stars (Table [l]). 

Johnstone k Guedel (2015) recently showed that the av¬ 
erage temperature of X-ray coronae (T cor ) is a weak function 
of X-ray flux F x : (T cor [MK]) = 0 . 11 (F x /[erg s^cm ” 2 ]) 0 "" 


(see also Tell eschi et al.|[2005]). Using their relation, the X- 
ray luminosities compiled in |Vidotto et ah] (2014b) and the 
stellar radii from Table [l] we estimate (T cor ) to be in the 
range between 2 and 3.6 MK for the stars in our sample. 
Naively, one could expect that the temperature at the base 
of the wind is related to the temperature of the closed X- 
ray corona (and this is the case for our Sun), but it is not 
clear if this relation is true for other stars. Therefore, in the 
absence of a stronger constraint, in our models, we adopt 
a wind base temperature of 2 MK, typical of stellar coro¬ 
nae of solar-type stars. We adopt 7 that is the same as the 


effective adiabatic index measured in the solar wind (Van 
Doorsselaere et ah||2011 ). For the base density, we adopted 
a value of lCPcm -3 . Ideally, observations of mass-loss rates 
of cool dwarf stars would allow us to place better constraints 
on the densities. However the lack of, or difficult-to-obtain, 
observational signatures of these winds make constraints of 
base density (or mass-loss rates) challenging to be obtained. 

To investigate how our results change with the change 
in the wind base density, we performed a stellar wind sim¬ 
ulation of HD 46375 that results in a mass-loss rate (M = 
2.9 x 10 -14 Mq yr -1 ) that is similar to the one observed in 


the solar wind (M = 2 x 10 “ 14 M® yr -1 )- Compared to the 
values of HD 46375 reported in Table [2] in this simulation, 
we found a mass-loss rate that is a factor of 6.5 smaller, J 
that is a factor 3 smaller and <f>o that is a factor 1.3 smaller. 

Locally, the hot-Jupiter HD 46375b experience a total 
external pressure whose average value (averaged over the 
longitude of the subplanetary point) is a factor of 5.6 smaller 
than the value presented in Table [ 3 ] Because vm is weakly 
dependent on p tot (tm oc p^/ 6 ), the value of vm we esti¬ 
mated before is smaller by a factor of only 1.3. In spite of 
the larger the cross-section of the planetary magnetosphere, 
the radio flux decreased by a factor of 2.3, caused by the 
decrease in the ram pressure [Eq. (|15|)]. 

Another parameter we have assumed in our models is 
the planetary magnetic field intensity. As discussed in Sec¬ 
tion 0 this is a quantity that has yet not been measured 
in exoplanets. Here, we adopted a magnetic field intensity 
which is similar to the value of Jupiter. We can also esti¬ 
mate how the magnetospheric size we presented in Fig. & 
would have changed if a different field strength were to be 
adopted. Because vm oc B^ 3 , a strength that is a factor of 
2 smaller would decrease the reported values of tm by 2 1 / 3 
(i.e., only 25%). In spite of that, this would not have signifi¬ 
cantly altered the computed radio flux (our radio flux model 
is weakly dependent on the planetary field strength; see dis¬ 
cussion in Vidotto et al. 2012), but would have decreased 


the frequency of the emission by a factor of 2 , making it not 
possible to be observed from the ground, due to the Earth’s 
ionospheric cut-off in frequencies. Indeed, one of the possibil¬ 
ities that exoplanetary radio emission has not been detected 
so far might be due to a frequency mismatch between the 


emission source and that of the search instruments (Bastian. 
Dulk k Leblanc|2000 ). 


6.2 Exoplanetary system conditions for 
detectability at radio wavelengths 

Because of the cyclotron nature of the magnetospheric radio 
emission, exoplanets with higher magnetic field strengths 
emit at higher frequencies, where the detection sensitivity is 
larger. For instance, an exoplanet with a magnetic field of 
about 40 - 50 G emits at the frequency range between 110 
and 140 MHz. The sensitivity of LOFAR at 100 to 200 MHz 


is roughly about 0.05 mJy (cf. Fig. 1 in GrieBmeier, Zarka 
k Girard|2011 ). This indicates that, except for HD 46375b, 


all the remaining exoplanets studied here could in principle 
be detectable with LOFAR, if their magnetic field strengths 
were about 40 - 50 G. Compared to Jupiter’ maximum field 
intensity, these field strengths are about ~ 3 times higher. 

We can also estimate what would be the required dis¬ 
sipated stellar wind power to generate detectable radio sig¬ 
natures from the exoplanets studied here. In this exercise, 
we take the same exoplanetary magnetic field assumed in 
Section [ 5 ] (i.e., B P:eq = 7 G). With such a magnetic field in¬ 
tensity, the frequency of emission is around 36 MHz, where 
the LOFAR sensitivity for a one-hour integration time is 
about a few mJy. The radio power calculated in Section [5] 
yielded values of about (1.6 - 5.6) X 10 25 erg s -1 . For a radio 
flux of a few mJy, the required radio power of the exoplan¬ 
ets studied here should be higher, on the range of (1.1 - 
2.1) x 10 27 erg s -1 . To have a radio power (or, equivalently, 
a wind kinetic power) that is roughly 2 orders of magni- 
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tude larger, the stellar wind characteristics need to change 
- either by increasing the density of the stellar wind or its 
velocity or both, as demonstrated next. 

From equations ( |15| ) and ( |16| ) , and assuming a ram 
pressure-dominated wind, one can show that 


p 2 Au 7 ■ 


8 Pi 


2 RIB* ’ 


(18) 


such that the ratio between the values required for a radio 
flux of about 3 mJy to the ratio of the values calculated at 
Section [5] is 

r [P 2 2 A AM J (3mJy) ~ f f fc(3mJy) V ~ (0.5 to 3.3) x 10 s . (19) 

[p Z Au r j(Sect. 5 ) \Bk(Sect. 5) J 

A very crude estimat^] then tells us that either the wind 
density needs to increase by a factor of at least ~ 300 to 


600 (i.e., the square root of the values derived in Eq. (19)) 
or the velocity requires an increase of a factor of at least 


~ 5 to 7 (i.e., the 7-th root of the values in Eq. (19)). Or, 
alternatively, density and velocity should both change such 


that they obey the relation (19) 


From Table |3| a 5 to 7 times increase in the wind ve¬ 
locity implies a relative velocity > 1200 km/s, which is 50% 
larger than the speed of the fast solar wind and 3 times larger 
than the slow solar wind speed. In terms of density, an in¬ 
crease of ~ 300 to 600 roughly implies a similar increase in 
mass loss rates and, from Table [ 2 ] this would result in M of 
at least (2.9 - 24) x 10 4 times the solar wind mass-loss rates. 
Such mass loss-rates are typical of very young stars, indicat¬ 
ing that exoplanets orbiting young Suns are more likely to 
produce detectable levels of radio fluxes (GrieBmeier et al. 
2005 Vidotto et al.||2010 ). 


7 SUMMARY AND CONCLUSIONS 

In this work we have investigated the interplanetary me¬ 
dia surrounding five hot-Jupiters, namely: HD 46375b, 
HD 73256b, HD 102195b, HD 130322b, HD 179949b. For 
that, we carried out 3D MHD stellar wind simulations, which 
incorporate as boundary conditions the surface magnetic 
held of the star reconstructed by Fares et al. (2012 2013) 
using the Zeeman Doppler Imaging technique. The global 
characteristics of our wind models are presented in Table [2] 
We then calculated the local characteristics of the stellar 
winds at the orbital radius of the hot-Jupiters, in order to 
characterise the interplanetary medium surrounding these 
exoplanets. In particular, we calculated the total pressure of 
the interplanetary medium and estimated what would be the 
size of planetary magnetospheres in case these hot-Jupiters 


2 Note that this approach is not a self-consistent one, because 
in this scenario the structure of the wind temperature, velocity 
and magnetic fields are not modified (i.e., we are assuming they 
remain unchanged as the structures derived in Section [4|. In the 
self-consistent approach, if either the density of the wind or its 
velocity are modified, one needs to solve the coupled MHD equa¬ 
tions to derive all the remaining quantities of the wind. However, 
this back-of-the-envelope calculation can give a rough estimate of 
how larger should the stellar wind power be in order for the radio 
emission to reach values above the sensitivity limit of a couple of 
mJy. 


had a magnetic held similar to Jupiter’s held. We found 
that magnetospheric sizes range between 4.1 and 5.6 R p and 
that they can vary by a few percent due to variations in the 
external environment of the planets, as they orbit around 
their parent stars. We also demonstrated that these plan¬ 
ets orbits are super fast magnetosonic, indicating that bow 
shocks should be formed around their magnetospheres. The 
bow shock orientations (i.e., the angle between the shock 
normal and the tangent of the circular orbit) are of interme¬ 
diate type, in which the shock forms at intermediate angles 
from the one of a shock facing the motion of the planet 
(‘ahead shock’) and the one connecting the star-planet cen¬ 
tres (‘dayside shock’). 

We also calculated the size of the auroral ovals of these 
planets. Inside these ovals, the planetary magnetic field lines 
are open, through which particles from the star and from the 
cosmos can penetrate as well as planetary atmospheric parti¬ 
cles can escape through polar flows. On average, the auroral 
ovals we calculated have a half-opening angle of about 25° 
to 29°, leaving exposed about 9% to 13% of the planetary 
area, which is a factor of ~ 2 larger than estimates for the 
Earth’s and Saturn’s auroral caps. Finally, we estimated the 
radio flux of these planets, using the analogy observed in the 
solar system, in which the radio emission from the magne¬ 
tised planets is correlated to the solar wind power. We found 
small radio fluxes ranging from 0.02 to 0.13 mJy, which 
should represent a challenge to be observed with present- 


day technology (e.g., LOFAR; GrieBmeier, Zarka V Girard 


2011 ), but could be detectable with higher sensitivity arrays, 


such as the SKA-low array system. Radio emission from sys¬ 
tems having closer hot-Jupiters, such as from r Boo b (ra¬ 
dio flux of the order of 0.5 - 0.9 mJy, Vidotto et al.||2012 ), 
HD 189733 b (0.5 - 1 mJy, calculated using the simulations 
from Llama et al. 2013| and the same model as presented 
here), or from nearby planetary systems orbiting young stars 
( | GrieBmeier et al.' 2005 Vidotto et al. 2010), are likely to 
have higher radio fluxes, presenting thus better prospects 
for detection of exoplanet ary radio emission. 
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